The 7*7 — * 7T° transition in a bound-state approach 
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We present a recent treatment of the 7r°7*7 transition form factor in the cou- 
pled Schwinger-Dyson and Bethe-Salpeter approach to stress the desirability of 
measuring it at CEBAF. 



The form factor T 7r o(— Q 2 ,0) for the transition ~f*(k)^(k') — > 7r°(p) (where 
k 2 = —Q 2 7^ is the momentum-squared of the spacelike off-shell photon 
7*), is usually thought to be almost completely understood on the basis of 
the perturbative QCD (pQCD) over the whole range of the momenta-squared 

covered by the new CLEO data gj (1.5 GeV 2 < Q 2 < 9 GeV 2 ). However, 
Radyushkin and Ruskov's Q QCD sum rule analysis indicates the presence 
of large nonperturbative contributions even for the largest of these presently 
accessible momenta. This motivated us to address in Ref. jjj the transition 
form factor in the coupled Schwinger-Dyson (SD) and Bethe-Salpeter (BS) 
approach in all spacelike-momentum regimes, from very large Q 2 all the way 
to Q 2 = 0. The Q 2 — limit corresponds to the tt° decay into two real 
photons (k 2 = k' 2 = 0) explained by the Abelian axial alias Adler-Bell-Jackiw 
(ABJ) anomaly, which is usually difficult to incorporate into a bound-state 
approach Fortunately, as shown already in Refs. |||| and used in the first 
analysis of T v o{— Q 2 ,0) in a SD approach jjj, such approaches which rely on 
the chirally well-behaved SD (and BS) equations for the light pseudoscalar 
mesons (tt, K,rj, ...) while respecting the Ward-Takahashi identities (WTI) of 
QED, reproduce (in the chiral and soft limit) the fundamental anomaly result 

2X0,0) = ^ (i) 

naturally and without any fine tuning. This is an advantage with respect to 
pQCD approaches, which have problems at low Q 2 . (E.g., see Ref. ||.) 



a The talk given by D. Klabucar at Semi-Exclusive Workshop, Jefferson Lab, May 19-22, 
1999. To appear in the proceedings of the Workshop. 



In the present coupled SD-BS approach the anomalous amplitude 

(|l|) is also obtained model-independently and without imposing any require- 
ments on the solutions for the ir° wave function or the quark propagator solu- 
tions (e.g., see This is exactly as it should be, since the anomaly - and 
its result ([l]) - must not depend on the internal structure of the pion. Sub- 
sequently, in Ref. ||, we showed that the SD-BS approach to modeling QCD 
provides a good description for both low and high values of Q 2 . We thereby 
obtained the behavior similar to the Brodsky-Lepage interpolation formula 
T^i-Q 2 , 0) = {1/4tt 2 U) I (1 + Q 2 /8ir 2 f 2 ) proposed 0] as a desirable behav- 
ior for T n o (— Q 2 , 0) as it reduces to the ABJ anomaly amplitude ([!]) at Q 2 = 0, 
while agreeing with the following type of leading behavior for large Q 2 : 

T>(-Q 2 ,0) = J (J = constant for large Q 2 ), (2) 

which is favored both experimentally M and theoretically (e.g., p|, [TT|.[T^] ) . 

In the coupled SD-BS approach, the BS equation for the pion bound- 
state qq vertex T^o (q, p) employs the dynamically dressed quark propagator 
S(k) — [A(k 2 )l/i — B(k 2 )]^ 1 , obtained by solving its SD equation. Then, in the 
case of light pseudoscalars, the qq bound states are simultaneously also the 
(pseudo-)Goldstone bosons of dynamical chiral symmetry breaking (D^SB). 

Following Jain and Munczek (T^,[l4|], we adopt the ladder-type approxima- 
tion employing bare quark-gluon-quark vertices but dressed propagators. For 
the gluon propagator we use an effective, (partially) modeled one in Landau- 
gauge @|3, given by G(~1 2 )(q^ - IH V jl 2 ) . (This Ansatz is often called 



the "Abelian approximation" |15|.) The effective propagator function G is the 
sum of the perturbative contribution Guv and the nonperturbative contribu- 
tion G IR : G(Q 2 ) = G UV (Q 2 ) + G IR (Q 2 ) , (Q 2 = -I 2 ) . The perturbative 
part Guv is required to reproduce correctly the ultraviolet (UV) asymptotic 
behavior that unambiguously follows from QCD in its high-energy regime. 
Therefore, Guv must be given by a s (Q 2 )/Q 2 , where the running coupling 
constant a s (Q 2 ) is well-known from pQCD, so that Guv is not modeled. As 
in Refs. ^,||,|l0), we follow Refs. and employ the two-loop asymptotic 

expression for a s (Q 2 ). For the modeled, IR part, we adopt from Ref. juj 
Gir(Q 2 ) = (167r 2 /3)aQ 2 e-^ 2 and a = (0.387 GeV)- 4 , fx = (0.510 GeV)- 2 . 

Our calculational procedures are detailed in our Refs. P,[^,[ic[| . We repro- 
duce the solutions of Ref. Jlij for the dressed propagator S(q), i.e., A(q 2 ) 
and B(q 2 ), as well as the solutions for the four functions comprising the 
pion bound-state vertex T^o. Actually, Ref. [14| employs the BS amplitude 
X 7r o(g,p) = S(q + p/2)T 7T a(q,p)S(q — p/2), which is completely equivalent. 

We assume that the 7r°7*7 transition (and other pseudoscalar meson — > 
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Figure 1: The GIA diagram for P° — > "y(*J»y(*J processes (P° = 7T°, 77, Tj' , T) c> Vb]- 



two photon transitions) proceeds through the pseudoscalar- vector- vector trian- 
gle graph (Fig. 1), and that we calculate the pertinent amplitude T^(fc, k') — 
£ Q ^fc Q fc^T w o(fc 2 ,fc' 2 ) as in Refs. using the framework 'advocated 

by, e.g., Refs. [fj-Q] in the context of electromagnetic interactions of BS bound 
states, and often called {e.g., by Ref. 0) the generalized impulse approxi- 
mation (GIA). We therefore use the dressed quark propagator S(q) and the 
pseudoscalar BS bound-state vertex Tp(q,p), as well as consistently dressed 
electromagnetic vertex T^lq', q), which satisfies the vector Ward-Takahashi 
identity (WTI) (q' - q) fi T' i (q', q) = S' 1 ^) - S~ l (q) . Namely, assuming that 
photons couple to quarks through the bare vertex 7^ would be inconsistent with 
our dressed quark propagator S(q), which contains the momentum-dependent 
functions A(q 2 ) and B(q 2 ). The bare vertex 7^ obviously violates the vector 
WTI, implying the nonconservation of the electromagnetic vector current and 
charge. Solving the pertinent SD equation for the dressed quark-photon-quark 
{ill) vertex T M is still an unsolved problem, and using the realistic Ansatze 
for T M still remains the only practical way to satisfy the WTI. The simplest 
particular solution of the vector WTI is the Ball-Chiu (BC) jl6j vertex 

^^( g ^ g ) = A4 g '^ 9 2 )^ + ^^{A_( g ^ g 2 )^^-i?_( ( z'^ g 2 )}, & 

where H±(q /2 ,q 2 ) = [H(q' 2 ) ± H(q 2 )], for H = A or B. It does not introduce 
any new parameters as it is completely determined by the dressed quark prop- 
agator S(q). In the SD-BS approach, this minimal WTI-satisfying Ansatz (||) 
is still the most widely used qq^f vertex (e.g., Refs. |p|,[^,p|,[lo|). A general WTI- 



3 



satisfying vertex can be written [jl6| as T M = T^ c + AT 11 , where the addition 
Ar M doesn't contribute to the WTI, since it is transverse, (<?' — q) f _ l AT fJ -(q' , q) = 
0. That is, AT fJ, (q',q) lies in the hyperplane spanned by the vectors 7~ M (q', g) 
(i = 1, ...,8) transverse to the photon momentum k — q' — q. Curtis and Pen- 
nington (CP) H^] advocated an Ansatz for AT^(q' 1 q) exclusively along the 
i = 6 basis vector: AW(q',q) = Tg{q', q)[A-(q' 2 , q 2 )/2d{q', q)]. The coeffi- 
cient of T§{q', q) can be chosen to ensure renormalizability |17|j , e.g., 



d±(q',q) = 



(q' 2 ±q 2 ) 2 + M 2 (q' 2 ) + M 2 (q 2 ) 



(4) 



where M(q 2 ) = B(q 2 ) / A(q 2 ) is the D^SB-generated dynamical mass function 
which in our case has the large-g 2 dependence fl in agreement with pQCD. 

The choice d = d_ corresponds to the CP vertex Ansatz T^ p 
use it in analytic calculations of X^o (— Q 2 , 0), which are possible for Q 2 



We 
= 



and Q — > oo. However, in the numerical calculations, which arc necessary for 
finite values of Q 2 ^ 0, we prefer the modified CP (mCP) vertex, r^ C p @, 
resulting from the choice d = d+, since it is easier to deal with numerically. 
Unlike the BC one, the mCP vertex is consistent with renormalizability just 
like the CP one. In the present context, the important qualitative difference 
between the BC vertex, and the CP as well as the modified, mCP vertex, is 



7 M when both q' 2 ,q 2 — ► ±oo, whereas T^ p (q\q) 



r 



and 



7 M as soon as one of the squared momenta tends to infinity. 



that r» c (q',q) 

T mCP 

Assuming the isospin symmetry, \ = X™ = XdS, IS the BS amplitude for 
both uu and dd pseudoscalar bound states, and Xn°(q,P) = x{q,p) A 3 /-\/2. 
Then, for 7r°, CIA yields (see Fig. I) the tensor amplitude 



q 



tr{r*(«-|,A- 



q - P -)S{k + q- P -) 



x T v {k + q-^,q+ \)x{q,P)} + (k~k',[x~ u), (5) 

leading to the analytical reproduction of the tt° — > 77 amplitude ([!]) with the 
CP and mCP vertices, in the same way as with the BC vertices. 

Let us stress that in the regime of asymptotically large Q 2 , we obtained || 
analytical results for the transition form factor in agreement with the behav- 
ior (||). Since the pion is light, k ■ k' w Q 2 /2 for large negative k 2 = —Q 2 
and k' 2 = 0. Taking into account the behavior of the propagator func- 



tions A(q 2 ) and B(q 2 ) (e.g. 



^]) we can then in Eq. Q approximate 
those quark propagators that depend on the photon momenta k and k', by 
their asymptotic forms: S(q - p/2 + k) w S(k - p/2) w -(2/Q 2 ){ft - rf/2) 
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and analogously for the propagator where k is replaced by k' . Although the 
relative loop momentum q can be large in the course of integration, its ne- 
glecting is justified because the BS amplitude x(l;P) decays quickly and thus 
strongly damps the integrand for the large q's. Ref. [|| then showed that 
T^{k,k') = ~(2/N c Q 2 )e^(k - k') x p a U, finally leading to Eq. §) with 
the coefficient J = 4/3. We thus found model independently that the asymp- 
totic behavior in the present approach (with qq^ vertices such as the bare 7 M , 
mCP and CP), agrees exactly with the leading term predicted by OPE ]T^ |. 

The asymptotic behavior oc 1/Q 2 obtained in the present approach is espe- 
cially satisfying when compared with the one known to result from the simple 
constituent quark model (with the constant light-quark mass parameter m u ), 
where T 7r o(-Q 2 ,0) oc (m 2 /Q 2 ) ln 2 (Q 2 /m 2 ) as Q 2 oo, which overshoots the 
data considerably because of the additional ln 2 (<3 2 )-dependence. 

It is important to note that this asymptotic behavior, T 7r o(— Q 2 ,0) ~ 
(4 / '3) /„■ / Q 2 , was obtained for the qq^ vertices which reduce to the bare one, 
7 M , as soon as one of the squared momenta tends to infinity, such as the renor- 
malizable CP and mCP vertices. Subsequently, Ref. |l8j showed that this 
asymptotic behavior must in fact hold for any renormalizable qq^f vertex. For 
the simplest dressed WTI-preserving vertex, the BC one, which reduces to the 
bare one only when momenta are large in both of its fermion legs, and which 
is not consistent with multiplicative renormalizability |l7j , the asymptotic be- 
havior is (4/3)/ T /Q 2 , where the quantity is given by the same Mandelstam- 
formalism expression as /„■ except that its integrand is modified by the factor 
[1 + A{q 2 )} 2 /4: |j. For our model solutions |j§, f n w 1.334/ T , so that the usage 
of the BC vertices causes the increase of J in Eq. |]irom J = 4/3 to J w 1.78. 

Of course, the model-independent asymptotic coefficient 4/^/3 is the one 
having the more fundamental meaning, resulting from the renormalizable qq"f 
vertices such as the CP or mCP ones, which have properties closer to the 
true vertex solution. Also indicative is the asymptotics found p] for such qqj 
vertices when both photons are off-shell, k 2 = —Q 2 << and k = —Q' 2 < 0: 

The usage of the BC vertex (^) would again modify this result by the sub- 
stitution /tt — > /tt. However, Eq. (^|) agrees with the leading term of the 
OPE result of Novikov et al. jl9| for the special case Q 2 = Q' 2 . Also, the 
distribution-amplitude-dependence of the pQCD approach cancels out for that 
symmetric case, so that T w o(— Q 2 , — Q' 2 ) in this approach {e.g., see |2Cj]), in 
the limit Q 2 = Q' 2 — > oo, exactly agrees with both our Eq. (||) and Ref. |fi1| . 
For that symmetric case, we should thus have even the precise agreement of 
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the coefficients irrespective of the description of the pion internal structure en- 
coded in the distribution amplitude. Obviously, this favors the renormalizable 
qq^ vertices, such as CP and mCP ones, over the BC vertex. However, the 
BC vertex may anyway be the one which is more accurate not only for the 
presently accessible Q 2 , but also for larger values before starting to fail. 




Figure 2: Our numerically obtained Q 2 T n o(— Q 2 ,0) and the data. 



For finite Q 2 , the transition form factor is evaluated numerically. Fig. 2 
compares so obtained Q 2 T n o (-Q 2 , 0) with the CLEO § and CELLO |lj data. 
The solid curve is obtained with the BC vertices, and the dashed curve with 
the mCP ones. The horizontal dashed line denotes 4/^/3, the asymptotic value 
for renormalizable qq"f vertices such as the mCP one. It barely touches the 
experimental error bars from below. Note that also Manohar Jl2] warned that 
his OPE approach (also yielding indicates the possibility of large correc- 

tions to his leading term, but possibly also to the pQCD asymptotic coefficient 
2/^, which is however exact in the strict ln(Q 2 ) — > oo limit |T[ . We get it lower 
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by the factor 2/3 because our present approach does not incorporate the ef- 
fects of the pQCD evolution. However, as Ref. || seems to indicate that effects 
other than the pQCD evolution (including the presently interesting dynamical 
dressing of quarks) may still play important (and maybe even dominant) role 
even at Q 2 -values larger than the presently accessible ones, one should consider 
seriously also nonperturbative approaches such as ours. On the other hand, 
high-precision measurements of T T o(— Q 2 , 0) can test (and give a hint on how 
to improve) the SD and BS model solutions which have so far been successful 
in fitting the low-energy hadron properties such as the meson spectrum. For 
example, our model choice jl4| somewhat overshoots (for both BC and mCP 

vertices) the present I>(-Q 2 ,0) data in tne region Q 2 < 4 GeV 2 . 

Once one has a solution for xiliP) that leads to the correct value of f v , 
the transition form factor is most sensitive on A(q 2 ), or, more precisely, on its 
values at small and intermediate momenta — q 2 , where A(q 2 ) is still appreciably 
different from 1 0,p2[. To illustrate what happens when the A(q 2 )-profile is 
decreased, let us enforce by hand the extreme, artificial case A(q 2 ) = 1. (To 
avoid confusion, we stress it is for illustrative purposes only, as we cannot 
have such a SD-solution in the adopted approach of Refs. |l^,|l^].) This leads 
to the curve traced on Fig. 2 by small crosses, pertaining to the usage of 
both BC and mCP vertices (as well as the CP ones), since A{q 2 ) = 1 makes 
^mCP ~ * ^bc- This cur ve reveals how the heights of the curves depicting 
Q 2 T 7r o (— Q 2 , 0) depend on how much the A(<7 2 )-profile exceeds 1. Obviously, for 
both the solid curve and the dashed one, the agreement with experiment would 
be improved by lowering them somewhat (at least in the momentum region 

Q 2 < 4 GeV), which could be achieved by modifying the model juj and/or its 
parameters so that such a new solution for A(q 2 ) is somewhat lowered towards 
its asymptotic value A(q 2 — > oo) — > 1. (Of course, in order to be significant, 
this must not be a specialized re- fitting aimed only at A(q 2 ). Lowering of A(q 2 ) 
should be a result of a broad fit to many meson properties, comparable to the 
original fit jlij. This, however, is beyond the scope of this paper.) 

By the same token, high precision measurements of T 7r o(— Q 2 ,0) can be 
especially helpful in obtaining information on which solutions for A(q 2 ) are em- 
pirically acceptable and which are not. Of course, measurements of T^o (— Q 2 , 0) 
give information on the integrated strength of A(q 2 ) rather than on A(q 2 ) it- 
self. However, since it is known that the form of that function must be a rather 
smooth transition (e.g., see [|J) from A(q 2 ) > 1 for q 2 near 0, to A(q 2 ) — > 1 in 
the q^-domain where QCD is perturbative, such measurements p3| would give 
a useful hint even about A(q 2 ) itself - namely about what solutions for A(q 2 ) 
one may have in sensible descriptions of dynamically dressed quarks and their 
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bound states. Therefore, the intermediate- momentum (Q 2 < 4 GeV 2 ) mea- 
surements of T^o (-Q 2 , 0) at Jefferson Lab, such as those proposed in Ref. Q , 
would be very desirable. Hopefully, the study of virtual Compton scattering 
at CEBAF can produce the data on T^o(~Q 2 , 0) as envisioned in Ref. P3[ . 
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